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The hierarchical scenario of structure formation describes how objects like galaxies and 
galaxy clusters are formed by mergers of small objects. In this scenario, mergers of 
galaxies can lead to the formation of massive black hole (MBH) binary systems. On the 
other hand, the merger of two MBH could produce a gravitational wave signal detectable, 
in principle, by the Laser Interferometer Space Antenna (LISA). In the present work, 
we use the Press?Schechter formalism, and its extension, to describe the merger rate of 
haloes which contain massive black holes. Here, we do not study the gravitational wave 
emission of these systems. However, we present an initial study to determine the number 
of systems formed via mergers that could permit, in a future extension of this work, the 
calculation of the signature in gravitational waves of these systems. 
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1. Introduction 

Recently, observational evidence for the existence of massive black holes in galaxies 
has been reported in the literature. On the other hand, the hierarchical scenario 
of structure formation describes how objects like galaxies and galaxy clusters are 
formed in the early universe. In this way, we can suppose that mergers of galaxies 
can lead to the formation of the massive black holes (MBH) observed in galaxies and 
their binary systems. Thus, the main goal of this work is to describe a method to 
determine the evolution of binary systems of massive black holes in the hierarchical 
scenario. To do that, we present in Sec. 2 a short review on the Press?Schechter 
formalism. In Sec. 3 we describe how to obtain the relation between the central black 
hole and the mass of the host dark halo. We also present the way to calculate the 
number of massive binary systems in Sec. 3. In Sec. 4 we present the main results, 
and finally in Sec. 5 we present our conclusions. Our models are obtained using the 
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following set of cosmological parameters: f2 m = 0.24, Slf, = 0.04, Ha — 0.76 and 
h = 0.73. 

2. Hierarchical scenario of structure formation: Press- Schechter 
formalism and its extension 

The hierarchical scenario of structure formation: the core of Press- SchechterlSl(P-S). 
formalism is that a dark matter halo leaves the linear regime when the mean density 
within a given volume is larger than a threshold level S c . In particular, Lacey and 
Cole ^proposed an extension of the P-S formalism based on the Brownian random 
wake of Bond et al. 1 ^ . The goal of this extension was to take into account the 
probability that a dark matter halo (henceforth halo), with mass Mi, has to merger 
with another halo with mass M2, for any redshift z, in order to form a new halo 
with mass M/ = Mi +M2. Fakhouri and MaP showed that, using the P-S formalism 
and its extension, the merger rate of haloes is given by: 

B(Mi,M f ,z) _ fl 1 
f(M fl z;P~S) V ™ 2 (Mi) 

where /(M, z;P — S) is the P-S mass function of dark haloes, a(M) is the variance 
and B(M\,Mf, z) is the merger of dark halos. 

3. Binary systems of massive black holes 

Wythe & Loeb0J proposed a model for the relation between central black hole and 
the mass of host dark halo. The authors consider that the central black hole (CBH) 
stop growing when the accretion reaches the Eddington luminosity. In particular, 
they consider that the circular velocity is equal the virial velocity, in this case, the 
mass of dark halo Mh, as a function of the CBH is (see Refs|31 [T]): 

where f2^ is the dark matter parameter at present time, £0 = 10 -5,7 and A c is the 
linear overdensity by virialization of a spherical pertubation il top-hat r '-like, that for 
ACDM isPl 

A c = 18tt 2 + 82[Q m (z) - 1] - 39[Q m (z) - l] 2 . (3) 

If we consider that the fraction e\ of dark haloes, at z < 10, having a central 
MBH,then we can obtain, using equations ([2]) and ([TJ, the following equation to the 
formation rate of massive binary systems: 

R(M bh ,i,Mbh,2,z) = eie2f(a)B(M h (M B n,i),M h (MBH,2),z), (4) 

where Mbh.i (i = 1, 2) is the mass of CBH, a is the separation of binary system at 
z and /(a) is the separation distribution function. 



efcr(Mi) 



dln{M x 



d5 c 
dz 



a 2 (Mi) 



-0/ * 



(1) 
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We considered in Eq. (0| that the central MBH quickly forms binary systems 
and which the fraction ti of these systems achieve the gravitational wave regime. 
On the other hand, for the separation distribution function we have: ^ 

'a\3/4 /a\3/2l 



f(a)da = ^ 



^. (5) 

a 



In the above equation we consider that the MBHs are formed in galaxy clusters. 
Thus, x represents the maximum separation of the components, and x is the typical 
dimension of galaxies clusters. 

The density number (ubh) of black hole binary systems obeys the conservation 
equation: 



dn 



BH 



dz 
~dl 



= (6) 



In a gravitational wave regime, the variation of separation with time is: 

dt oc a 
Finally, the number of systems, as a function of z and observed frequency, v b s , 

is: 

dV da 
dz dv ohs ' 



N sys = -n BH — x— • (8) 



with: 



[G(M BHS + M BHa )] 1/2 a- 5 / 2 (l + z) (9) 



dv obs 2tt 
where dV/dz is the comovel volume. 

4. Numerical results 

In this work, we use the Lax-Wendroff schema^ to obtain the numerical solution 
of equation ([6]). Figure [1] on the left hand side, shows the number of systems, by 
frequency by redshift, of binary system formed by black holes of mass M B h.i = 
10 5 Mq and Mbh,2 — 0-WMbh,i and, on the right hand side, the total number of 
systems into the mass range 10 4 M Q < Mbh,i < 1O 7 M and 0.1M B h,i < M B h.2 < 
M B h,i for different values of ei and £2. In both cases, we assumed that the separation 
range is 3(r s / li i + r s h.2) < a < 100(r s ^ i + ^2), with r s h.i is the Schwardschild 
ratio of black hole i and we assumed x — 1.5 Mpc. Note that the function Nbh 
peaks at z w 2. 

5. Conclusion 

We presented a different method to calculate the number of binary systems of mas- 
sive black holes using the Press-Schechter formalism. It is important to emphasize 
that the total number of systems obtained here takes into account all systems. This 
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Fig. 1. The left hand side is showed the number of systems as a function of observed frequency 
and z with ei = £2 = 1 and black holes with mass Mbh,1 = 10 s and Mbh,2 = O.XMbh,1- ln the 
right hand right, is presented the total number of binary sistems of massive black holes into the 
mass range 1O 4 M < M B h,i < 10 7 M Q and 0.1M B h,i < MbH,2 < MbH,i- On the top it was 
assumed ei = 62 = 1 and on the botton it was considered ei = 0.01, £2 = 0.001, e\ = 0.09 and 
£2 = 0.0001. 

is a different result when compared, for example, for the work of Wythe & Loeb 
- who obtained their results only for systems within the LISA range. The numeri- 
cal method used here is stable and it produces a smooth function for the density of 
binary systems. The same method was used by Banerjee & Ghosh© or the calculus 
of the formation of binary systems in globular clusters of stars. 
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